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ACTIVE CONTROL OF A FLEXIBLE BEAM 


Because of inherent low damping and high flexibility, large space structures 
may require some form of active control of their dynamics. Because of the apparent 
inability to accurately model the dynamics of these structures, methods for 
parameter adaptive control are now being developed at Langley. The process 
currently being studied is shown in the block diagram below. This approach uses a 
digital computer to process discrete sensor data, identify modal parameters, 
calculate modal control gains, and then convert the modal forces to real forces. 

The last two blocks are the topic of this presentation. Some of the problems con- 
sidered are: Cl) the possibility that there may be many modes to control with 

limited amounts of hardware, and (2) the required accuracy of identified structural 
parameters. 


BACKGROUND: 

- NEED TO CONTROL FLEXIBLE MOTION OF LARGE SPACE STRUCTURES 

- ABILITY TO ACCURATELY MODEL THE DYNAMICS OF THESE STRUCTURES IS UNCERTAIN 

- THEORY NOW BEING DEVELOPED FOR PARAMETER ADAPTIVE CONTROL OF THESE STRUCTURES 



PROBLEMS: 

- POSSIBLY MANY MODES TO CONTROL 

- LIMITED HARDWARE (COMPUTATION, SENSORS, ACTUATORS) 

- HOW WELL STRUCTURAL PARAMETERS MUST BE IDENTIFIED FOR CONTROL 
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RESEARCH TASK 


The specific research task was to design a digital control scheme to suppress 
vibration of a homogeneous free-free beam. A digital computer simulation algorithm 
was then used to test (1) the effects of controlling more modes than available 
actuators, and (2) the sensitivity to identified structural parameters. 


• DESIGN A DIGITAL CONTROL SCHEME TO SUPPRESS VIBRATION OF A 
HOMOGENEOUS FREE-FREE BEAM 

• EXAMINE EFFECT ON STABILITY OF: 

- FEWER ACTUATORS THAN CONTROLLED MODES 

- ERRORS IN STRUCTURAL MODEL PARAMETERS 


t TEST WITH AN EXISTING SIMULATION ALGORITHM 



MATHEMATICAL MODEL 


The continuous beam was modeled by using the SPAR finite element algorithm 
which generates mode shapes and frequencies. These were used to write a modal 
representation of the beam dynamics which was used to design the control gains. 




t 


CONTINUOUS MODEL 



FINITE ELEMENT MODEL OF BEAM 


MODE SHAPES AND FREQUENCIES FOR 25 ELEMENT MODEL ARE 
GENERATED BY SPAR 
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t MODAL REPRESENTATION 

x = [E] q 
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- SET OF 50 UNCOUPLED 2ND-0RDER SYSTEMS 


• CONTROL SYSTEM DESIGNED USING THIS MODAL REPRESENTATION OF THE 
STRUCTURAL DYNAMICS 
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DISCRETE TIME MODEL 

In order to simplify the digital simulation and prepare for eventual digital 
implementation, the modal equations of motion were discretized. This results in 
the scalar equation which shows the present modal amplitude to be a function of the 
past two amplitudes and the past two controls. There is a discrete time transforma- 
tion analogous to the Laplace transform which results in a characteristic polynomial 
in z. The roots of this polynomial may be plotted in the complex plane with 
stability represented by magnitudes of less than 1. Analysis of the control system 
will be done primarily in this z-plane. 


• NECESSARY FOR: 

- DESIGN OF DIGITAL SIMULATION 

- EVENTUAL DIGITAL IMPLEMENTATION 

• DISCRETE EQUATION OF MOTION: 

q (k) = A l q (k-l) + A 2 q ( k-2 ) + B i u ( k- 1 ) + 8 2 u (k-2)' 
where A = 1 : (w ,t), t = sampling interval, 

B = f(w ,e,T), e = mode shapes. 



(k-2)t ( k - 1 ) t kr 


• DISCRETE TIME TRANSFORM (ANALOGOUS TO LAPLACE TRANS-) RESULTS Time 
IN A CHARACTERISTIC POLYNOMIAL IN Z: n 
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CONTROL DESIGN APPROACH 


The modal control design approach is to choose desired closed- loop roots from 
which the modal controller gains can be calculated. The modal control forces may 
be calculated directly, and the actual control forces can be calculated using a 
pseudo-inverse . 


t CHOOSE DESIRED CLOSED LOOP ROOTS FOR EACH MODE 

• CALCULATE MODAL CONTROLLER GAINS 

• CALCULATE MODAL CONTROL FORCES 

• CONVERT MODAL FORCES TO ACTUATOR FORCES 
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DIGITAL CONTROL OF ONE MODE 

The control of one mode is achieved by using the minimum order control law 
required for pole placement. This is of the same form as the modal equation of 
motion. The closed loop controller has a fourth-order characteristic equation as 
shown. The coefficients of this equation are determined from the desired closed- 
loop roots and are functions of the mode and control coefficients in the plant and 
control modal equations. The control objective is to achieve the desired roots 
defined by (a,b,c,d) by solving for the controller gains (C^, C^. D^, D^) . 


• MINIMUM ORDER CONTROL LAW FOR POLE PLACEMENT 

u (k) ’ C 1 q ( k- 1 ) + C 2 q (k-2) + D 1 U ( k - 1 ) + °2 U ( k-2 ) 

- NOTE SAME FORM AS PLANT MODEL 

[q (k) * A 1 q (k-l) + A 2 q ( k- 2 ) + B 1 u (k-l) + B 2 u (k-2) ] 

- C'S & D's ARE FOUR CONTROL GAINS 

• CLOSED LOOP CONTROLLER 



- HAS DISCRETE TIME CHARACTERISTIC EQUATION 

4 3 2 

z + az +bz +cz + d = 0 

where (a ,b ,c ,d } = f(A ,B ,C ,D ). 

• CONTROL OBJECTIVE: 

- ACHIEVE DESIRED CLOSED LOOP ROOTS (AS DEFINED BY a,b,c,d) 
BY CALCULATING THE CONTROL GAINS, 
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CLOSED-LOOP CHARACTERISTICS 

The control problem is shown in the block diagram below. The parameter 
uncertainties affect the calculation of the modal controller gains and a pseudo- 
inverse results in a least squares type error in the actual forces applied to the 
beam. A typical set of closed-loop roots is shown in each of the two complex 
plane plots. The design criterion is to place the roots of the controlled modes on 
a constant damping line in the s-plane. This line maps onto a constant- radius 
circle on the complex z-plane, with uncontrolled roots on the unit circle. 


o CONTROL PROBLEM 



o REPRESENTATION OF CLOSED LOOP DYNAMICS 
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RESULTS - CASE 1 


A baseline set of results is shown below. Exact parameters are used to 
calculate the control gains and five modes are controlled with five actuators. 
Note that all roots are calculated exactly so that all modal amplitudes have the 
same decay envelope. Also, mode six, which is not controlled, does show "minor" 
excitation and continues to "ring" after control to the other modes is stopped. 
This is not evident from the figure. 


* CASE I : EQUAL NUMBER OF ACTUATORS AND CONTROLLED MODES 
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RESULTS - CASE II 


This next caae illustrates an attempt to control more modes than available 
actuators. Here it is no longer possible to solve exactly for t he desired roots, 
and one root does become unstable, as shown by the increasing amplitude of the 
fourth mode . 


• CASE .1 1 : FEWER ACTUATORS THAN CONTROLLED MODES - UNSTABLE 


5 CONTROLLED 


Z-PLRNE 
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RESULTS - CASE III 


In the third case it is shown that attempting to control more modes than 
available actuators does not necessarily mean the system will be unstable. Note 
here that the calculated roots are closer to the desired roots and that, while the 
fourth mode is near the unit circle, it now is slightly stable. This can be 
explained by considering that the additional actuator and mode provide one more 
data point for the least squares fit. 


t CASE III: FEWER ACTUATORS THAN CONTROLLED MODES - STABLE 
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PARAMETER ERROR SENSITIVITY 


The results of the parameter error sensitivity study are shown below. Errors 
was placed on the modal frequency and damping parameters of the inodes selected for 
control. Experimental results show that the parameter estimator may have errors on 
the order of ±5 percent. The locus of roots calculated using a parameter error 
range of ±20 percent is plotted below. The magnitude of the error at which the 
inodes became unstable is summarized in the table. Note that as the mode number 
increases, the sensitivity decreases for the first three modes. The high sensi- 
tivity of mode six in case III is unexplained at this time. 


• ERROR ON MODAL FREQUENCY AND DAMPING PARAMETERS ONLY 

q (M = lIl q (k - i ) + dL q (k-2) * B 1 u {k-l) + u ( k - 2 ) 

• NO ERROR ON RIGID BODY PARAMETERS 

• BASED ON EXPERIMENTAL PARAMETER ESTIMATE ERRORS OF + 5Z 

5 CONTROLLED 
5 ACTUATORS 

Z-PLRNE 




Parameter error at instability 

CASE 

MODE 3 

MODE 4 ( MODE 5 

MODE 6 

I 5 CONT. 
5 ACT, 

-10 % 
+8 * 

-12% 

-14% 

X 

I I 5 CONT. 

R ACT. 

-6 % 
+ 14% 

0.0 % 

-14 % 


I I I 6 CONT. 
5 ACT. 

-8 % 
+8 % 

-12 % 

-16% 

-2 % 
+2 % 
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CONCLUSIONS 


1. LIMITED ACTUATORS: 

- THIS DESIGN PROCESS YIELDS UNDESIRABLE CLOSED LOOP DYNAMICS 
WHEN THE NUMBER OF CONTROLLED MODES EXCEEDS THE NUMBER OF 
AVAILABLE ACTUATORS, 

2. PARAMETER ERROR: 

- ERRORS WITHIN THE RANGE OF EXPERIMENTAL RESULTS CAN CAUSE 
INSTABILITY. 

- CONTROL SYSTEM MUST BE MADE MORE TOLERANT OF PARAMETER ERROR. 
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